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ABSTRACT 

We discuss the cooling history of the low-mass, helium core white dwarfs in short 
orbital period millisecond pulsars PSR J0751+1807 and PSR J1012+5307. The re- 
vised cooling age estimated by Alberts et al. agrees with the age estimation for PSR 
J1012+5307, removing the discrepancy between the spin-down age and the cooling 
age. However, if we accept this model then the helium white dwarf in the binary pul- 
sar system PSR J0751+1807 must be much hotter than is observed. We propose that 
this discrepancy may be resolved if, after detachment of the secondary star from its 
Roche lobe in PSR J0751+1807, the star loses its hydrogen envelope due to pulsar 
irradiation. When hydrogen burning stops, the white dwarf will cool down much more 
quickly than in the case of a thick hydrogen envelope with a hydrogen burning shell. 
We discuss several possibilities to explain different cooling histories of white dwarfs in 
both systems. 
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1 INTRODUCTION 

It is accepted that formation of a binary millisecond (or 
recycled) pulsar with a low-mass companion may be ex- 
plained as the end-point of close binary evolution in which 
an old pulsar is spun-up by accretion from the red giant 
(secondary) (see Bhattacharaya & van den Heuvel 1991, 
' for review). Once the secondary overflows its Roche lobe, 
it starts transferring mass and angular momentum to the 
neutron star. During the long term accretion phase the neu- 
tron star is generally assumed to be spun-up to the equilib- 
rium period (Peq)- The equilibrium period can be written as 
(Bhattacharaya & van den Heuvel 1991): 

where B 9 = B/10 9 G, R 6 = i?„ s /10 6 cm and M ns are the 
magnetic field strength, radius and mass of the neutron star 
respectively and MEdd is the Eddington mass accretion rate. 
After the whole envelope of the red giant has been trans- 
ferred, the binary consists of a low-mass star with a helium 
core (progenitor of a low-mass helium white dwarf) and a 
fast rotating neutron star (millisecond pulsar). At this point 
the pulsar spin period starts to change due to magneto- 



dipole radiation, with a characteristic time t elapsed from 
this epoch 

t = TX[ ^h ][1 - ( T r ~ 1] > (2) 

where t—P/(2P), Po and P are the initial and current spin 
periods, and P is the spin period derivative, n is the braking 
index (=3) and if Po "CP then t=r. It is possible to estimate 
the age t from the cooling age of the white dwarf. The latter 
can be calculated from the observed parameters of the (T e //, 
log g). This gives a unique opportunity for determination of 
neutron star age which is independently of the spin-down 
theory of pulsars. 

In this short note we shall discuss cooling properties 
of helium white dwarfs in short period binary systems with 
millisecond pulsars, PSR J0751+1807 and PSR J1012+5307. 



2 OBSERVATIONAL DATA FOR 
MILLISECOND PULSARS 

For both systems, the characteristic time t (eq (2)), to be 
compared with the cooling age of the white dwarf, is as- 
sumed to be equal to the spin-down age r of the neutron 
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star. This requires Po ~ Peq -C P. The last restriction will 
be relaxed in the discussion section (par. 6). 



2.1 PSR J0751 + 1807 

This system consists of a 3.5 ms pulsar in a 6 h orbit with a 
~ 0.15 Mq companion (Lundgren, Zepka & Cordes 1995). 
We estimate the magnetic field strength of neutron star from 
the relation 



B = 3x 10 19 (PP) 1/2 . 



(3) 



For the observed P and P (8 x 1CT 21 ss _1 ) we have B 
~ 1.7x10 G, and the spin-down age of the neutron star, r 
~ 8 Gyr. For this age and M wc i ~ 0.15 Mq a simple white 
dwarf cooling model by Mestel (1952) gives L ~ W~ 43 Lq. 
Lundgren at al. (1996) estimated a lower bound for the vi- 
sual magnitude, m v >23.5, which requires T e // < 9000 K 
and an age greater than 7.94 x 10 s years. From their data a 
main-sequence companion is ruled out down to an M5 star 
(M ~ 0.10 Mq), since such a star would have been detected 
at the distance of the pulsar. Therefore the pulsar compan- 
ion in this system is clearly a white dwarf. 



2.2 PSR J1012+5307 

PSR J1012+5307 consists of a millisecond pulsar (with P = 
5.256 ms) in a 14.52 h orbit with a ~ 0.15 Mq companion 
(M ns =lA Mq, i = 60°). P=1.46xl0 _20 ss _1 leads to the 
spin-down age, r ~ 7.0 Gyr (Nicastro et al. 1995, Lorimer 
et al. 1995), and B ~ 2.6 xlO 8 G. Van Kerwijk et al. (1996) 
and Callanan et al. (1998) estimated the mass of the he- 
lium white dwarf to be 0.16 ± 0.02 Mq. Taking the obser- 
vationaly determined gravity and effective temperature for 
the white dwarf, and adopting the cooling model of Iben & 
Tutukov (1986), the cooling age has been estimated as 0.3 
Gyr (Lorimer et al. 1995). The discrepancy between these 
two age estimations is solved if, in the white dwarf cooling 
models, a thick hydrogen envelope with stationary hydro- 
gen shell burning is taken into account (Alberts et al. 1996, 
Driebe et al. 1998), rising the cooling age to several Gyr. 



3 THE EVOLUTIONARY CODE 

The evolutionary sequences we have calculated consist of 
three main phases: 

(i) detached evolution lasting until the companion fills its 
Roche lobe, 

(ii) semi-detached evolution, 

(iii) a cooling phase of the white dwarf (the final phase dur- 
ing which a system with a ms pulsar + low-mass helium 
white dwarf is left behind). 



3.1 Detached phase 

The duration of the detached phase is somewhat uncertain; 
it may be determined either by the nuclear time-scale or by 
the much shorter time-scale of the orbital angular momen- 
tum loss owing to the magnetized stellar wind. 



3.2 Semi— detached phase 

In our calculations we assume that the semi-detached evo- 
lution of a binary system is non-conservative, i.e. the total 
mass and angular momentum of the system are not con- 
served. The formalism which we have adopted is described 
in Muslimov & Sarna (1993). We introduce the parameter, 
/i, characterizing the loss of mass from the binary system 
and defined by the relations, 



M = M 2 fi and Mi = -M 2 (l - /i), 



(4) 



where M is the mass-loss rate from the system, M2 
is the rate of mass-loss from the secondary star, and Mi 
is the accretion rate onto the neutron star (primary). The 
matter leaving the system will remove angular momentum 
according to the formula 



= h 



M 2 M 



yr 



(5) 



where M = M1+M2, Mi and M2 are the masses of the 
neutron star and secondary star respectively. Here we have 
introduced an additional parameter, f 2 , which describes the 
efficiency of the orbital angular momentum loss from the 
system due to a stellar wind (Tout & Hall 1991). In our 
calculations we assume /i=l and /2 = 1. 

We also assume that the secondary star, possessing a 
convective envelope, experiences magnetic braking (Mestel 
1968; Mestel & Spruit 1987; Muslimov & Sarna 1995), and 
as a consequence of this, the system loses its orbital angular 
momentum. For magnetic stellar wind we used the formula 
for the orbital angular momentum loss 



= -3 x 10" 



r M 2 Rl 

M1M20 5 



yr 



(6) 



where a and R 2 are the separation of the components 
and the radius of the secondary star in solar units. 

In the semi-detached phase we have also included the 
effect of illumination of the secondary star by the millisecond 
pulsar. In our calculations we assume that the atmosphere 
of the secondary star is heated by the hard (X-ray and 7— 
ray) radiation from the pulsar (Muslimov & Sarna 1993). 
The effective temperature, T e //, of the companion during 
the illumination stage is determined from the relation 



+ P M =4naR 2 2 T* ff , 



(7) 



where L in is the intrinsic luminosity corresponding to 
the radiation flux coming from the stellar interior, a is the 
Stefan-Boltzmann constant and R 2 is the radius of the sec- 
ondary. Put is the millisecond pulsar radiation that heats 
the photosphere, given by 



'ill = fz (y 



2a J 



(8) 



where L rot is "rotational luminosity" of the neutron star 
due to magneto-dipole radiation (plus a wind of relativistic 
particles) 

2tt^ 4 



Lrot — -^B 2 R r [ 



(9) 



in the above formula Rns is the neutron star radius, fz 
is a factor introduced to account for the efficiency of trans- 
formation of the irradiation flux into thermal energy (in our 
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case we take / 3 = 2 x 1(T 3 Muslimov & Sarna (1993)). Note 
that in our calculations the effect of irradiation is formally 
treated by means of modification of the outer boundary con- 
dition, according to the above relation. 



3.3 Cooling phase 

For systems in the white dwarf cooling phase, we take into 
account two processes: 

(i) the loss of orbital angular momentum due to the emission 
of gravitational radiation (Landau & Lifshitz 1971), 

(ii) an irradiation induced stellar wind, with corresponding 
loss of orbital angular momentum estimated as 



1 + q M 2 



(10) 



where q is the mass ratio (q = M\/M 2 ) and M wind is a 
mass loss of an irradiation induced stellar wind defined by 
the simple van den Heuvel & van Paradijs (1988) model 



Jh_ (R 2 \ 2 2R b ns a /2tt\ 
1 GM 2 \ 2a) 3c 3 V P ) 



(11) 



where / is efficiency factor (from to 1) and a is the 
separation between the stars in solar units. In our calcu- 
lations we assume three values of /: 0, 0.1 and 0.5 (very 
efficient irradiation). 



3.4 The code 

The models of the stars filling their Roche lobes were com- 
puted using a standard stellar evolution code based on 
the Henyey-type code of Paczyhski (1970), which has been 
adapted to low-mass stars. Convection is treated with the 
mixing-length algorithm proposed by Paczyhski (1969). We 
solve the problem of radiative transport by employing the 
opacity tables of Iglesias & Rogers (1996). Where the Igle- 
sias & Rogers (1996) tables are incomplete, we switch to 
opacity tables of Heubner et al. (1977). For temperatures 
less than 6000 K we use opacities calculated by Alexander 
& Ferguson (1994) and Alexander (private communication). 
Finally, we assume a Population I chemical composition for 
the secondary (X=0.7, Z=0.02). 



4 EVOLUTIONARY CONSIDERATIONS 

To produce short orbital period systems composed of a low- 
mass helium white dwarf and a millisecond pulsar, it is nec- 
essary that the low-mass secondary (red giant) fills its Roche 
lobe below the so called bifurcation period, and above the 
boundary period (Ergma & Sarna 1996, Ergma, Sarna & 
Antipova 1998). Ergma et al. (1998) and recent calculations 
by Sarna, Ergma & Antipova (2000) have shown that af- 
ter detachment of the secondary star from its Roche lobe, 
a rather thick hydrogen layer is left on the top of the he- 
lium core. The same result has been obtained by Driebe 
et al. (1998). This thick hydrogen layer will greatly influ- 
ence the cooling history of the very low-mass helium white 
dwarf (< 0.25 Mq). Alberts et al. (1996) were the first to 
argue that because white dwarfs with M md <Q.2Q Mq do 
not show thermal flashes (as found by Webbink, 1975), the 



amount of hydrogen in the envelopes must remain high, re- 
sulting in a much longer phase of significant hydrogen burn- 
ing. Thus cooling times for these objects are considerably 
extended. This removes the discrepancy between the cool- 
ing and spin-down ages for the millisecond binary systems 
PSR J1012+5307 (Alberts et al. 1996, Driebe et al. 1998, 
Sarna et al. 1999,2000). 

But now the question arises for the PSR J0751+1807. 
The evolutionary history for both pulsar systems are very 
similar, as argued by Ergma & Sarna (1996). This means 
that in PSR J0751+1807 the helium white dwarf must 
also be in the stable, long-lasting hydrogen-burning phase, 
which leads to long cooling times. However, observations 
show that the white dwarf in this system is undetectable, 
and really very old. 

To solve this discrepancy, we propose that irradia- 
tion of the secondary star by the millisecond pulsar PSR 
J0751+1807 causes mass loss after detachment from the 
Roche lobe. After detachment the star with a hydrogen 
burning shell evolves from the red giant region towards 
higher effective temperatures. If we accept the simple ir- 
radiation driven mass loss model proposed by van den 
Heuvel & van Paradijs (1988), then using eq. 11, for R 2 = 
0.56 R e (see point 2 on Fig. 1), a = 2 Rq, M 2 =0.15 Mq, 
Rns = 10 km and / = 0.1 (efficiency factor), we have 
M w i n d ~ 4xl0 -10 MQ-yr^ 1 . Since the evolutionary time- 
scale to evolve is rather long (at least several hundred mil- 
lion years), about 0.01 Mq of the hydrogen envelope will be 
lost. As a result, the hydrogen burning shell will switch off, 
and the white dwarf may cool down very quickly. 



5 SIMPLE COOLING MODEL FOR THE 
WHITE DWARF IN PSR J0751 + 1807 

To illustrate how the mass loss due to irradiation will in- 
fluence the cooling evolution of the pre-white dwarf (after 
detachment from Roche lobe), we computed the evolution- 
ary sequence for a system of 1.5 + 1.4 Mq, with P or b,i 
(RLOF)=1.07 d (RLOF - Roche lobe overflow) and Z=0.02 
(Fig. 1). After detachment from the Roche lobe, mass loss 
at a rate estimated from eq. 11 has been included (Fig. 2). 

The mass loss phase leads to a change of the orbital an- 
gular momentum, which leads to significant changes in the 
orbital period evolution (Fig.3). Figure 4 shows the evolu- 
tion of the visual magnitude (mv) of the progenitor of the 
white dwarf during the cooling phase, with wind mass loss 
included. 

In the H-R diagram (Fig. 1) the evolution of this sys- 
tem is presented with and without mass loss. The decrease 
of the mass of the hydrogen envelope leads to a significant 
decrease of the hydrogen shell burning efficiency. Therefore, 
if compared to the evolution of a white dwarf with stationary 
hydrogen shell burning, the cooling time is decreased. For 
example, when log L/L© = -2.0, the cooling age without 
mass loss is 8.4 x 10 9 yrs, with moderate mass loss (f=0.1) 
is 5.07 x 10 9 yrs, and with very efficient mass loss (f=0.5) is 
2.17 x 10 9 yrs. Due to angular momentum loss by the irradi- 
ation induced wind, the orbital period and separation of the 
system (after detachment of the secondary from the Roche 
lobe) initially increases (Fig. 3) . Because of this increase in 
the separation, and shrinking of the white dwarf (see eq. 11) 
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4.2 4 3.8 3.6 

log T eff 

Figure 1. HR diagram with evolutionary tracks of a 1.5 Mq 
star with 1.4 Mq neutron star. Different phases of evolution 
are: detached evolution before filling the Roche lobe (0-1), semi- 
detached evolution (1-2), and detached evolution during the cool- 
ing of a white dwarf (2—3). Phase (2—3) is plotted for three differ- 
ent cases; normal cooling (thick solid line), cooling with irradia- 
tion induced wind (f=0.1 - thick dashed line) and very efficient 
irradiation induced wind (f=0.5 - thick dash-dotted line). The 
thin solid line shows the ZAMS as calculated using our code. 
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5 6 7 8 9 10 



log t cool 

Figure 3. Evolution of the orbital period during the cooling 
phase for the three cases considered. The meaning of the lines 
is the same as for Fig. 1. The thin dashed line shows location, 
at which the surface effective temperature of the white dwarf is 
equal to 9000 K. 
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Figure 2. The dependence of the mass loss rate as a function of 
time in two different phases of evolution: semi-detached evolution 
(phase (1—2) — solid line) and cooling with irradiation induced 
wind. For the wind phase we show two different efficiencies of the 
wind: f=0.1 (dashed line) and f=0.5 (dash-dotted line). 



the mass loss rate (Fig. 2) decreases by two orders of magni- 
tude, making this mechanism of orbital angular momentum 
loss insignificant. The orbital period starts to decay due to 
emission of gravitational waves. 

Figure 4 gives the limit for acceptable my from theo- 
retical cooling tracks. We can exclude normal cooling (with 
shell hydrogen burning) due to a too low my. Two other 
cases of cooling with irradiation induced wind give the right 
range of visual magnitude. For the upper limit of pulsar age 
(spin-down age ~ 8 Gyrs) my and T ef f are (22.4, 7200K) 
for / = 0.1 and (24.7, 4700K) for / = 0.5. 

One can ask why we do not observe circumbinary mat- 
ter inside PSR J0751+1807, as we do in PSR 1957+20 or 
PSR 1744-24A (Lyne 1995). We propose that after having 
lost most of its hydrogen envelope (during the high mass loss 
phase) ,the star hydrogen burning stops and relaxes to a low- 
mass helium white dwarf. Irradiation of the white dwarf by 
the pulsar may be significant, but it does not lead to change 
in the internal structure of the envelope (as it does for a 
low-mass star with a convective envelope and shell nuclear 
burning). Therefore, no hot extended envelope is created 
and no mass loss induced. 
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Figure 4. Evolution of the visual magnitude during the cooling 
phase for the three cases considered. The meaning of lines is the 
same as for Fig. 1. The thin horizontal dashed line shows the 
Lundgren et al. (1996) lower bound for visual magnitude. The 
two thin vertical dashed lines show lower and upper limits for the 
millisecond pulsar age. 

6 DISCUSSION 

The first question is what is the real value of the initial 
pulsar period Po for both systems? Is it equal to P eq ? To 
explain the discrepancy between the cooling and spin-down 
ages of PSR J1012+5307 Burderi et al. (1996) have accepted 
the so called accretion-induced magnetic field decay model 
(see Bhattacharya & Srinivasan 1995 for review). In this 
case, at the end of the accretion phase, field decay takes place 
on a shorter time-scale than spin-up, and spin equilibrium 
is no longer achieved. 

The second point is connected with the discussion of 
Webbink's (1975) results, who found that there is a hydro- 
gen burning in the thick envelope (of a few x10~ 3 Mq) of 
low-mass helium white dwarfs (< 0.2 Mq). Recently three 
different groups, Albert et al. (1996, using Eggleton evolu- 
tionary program), Driebe et al. 1998 (Kippenhahn program) 
and we (Sarna et al. 1999,2000, Paczyriski program) have 
obtained the same results: small mass helium white dwarfs 
(< O.2M0), after detachment from the Roche lobe, have a 
massive hydrogen layer (a few xl0~ 2 Mq) in which station- 
ary hydrogen burning continues. This keeps these low-mass 
helium white dwarfs hot and prolongs their cooling time. 

However, there is a difference between Alberts et al. 
and our+Driebe's results, since we and Driebe found hydro- 
gen flashes for M> 0.2 Mq and Alberts et al. did not. We 



show that the transition from a cooling white dwarf with 
a temporarily stable hydrogen-burning shell to a cooling 
white dwarf in which almost all residual hydrogen is lost in 
a few thermal flashes (via Roche-lobe overflow) occurs be- 
tween 0.183-0.213 Mq, which well agrees with the cooling 
scenario proposed by Burderi et al. (1998). 

Therefore we may accept that PSR J1012+5307 is in 
a long term hydrogen shell burning stage and this may be 
the reason why the white dwarf is too hot (accepting that 
spin-down age of 7 Gyr, which value could be even larger 
if the pulsar has a significant transverse velocity (Hansen 
& Phinney 1998), is almost equal to the cooling age). For 
PSR J0751+1807 the irradiation driven mass loss reduces 
the hydrogen envelope mass, the hydrogen burning deceases 
and the white dwarf is quickly cooled. 

Taking into account the effect of illumination, P eq is 
- 1.6 x 10~ 3 s (Ergma & Sarna 1996). Therefore, if the pul- 
sar in PSR J1012+5307 has been spun up to P eq , it also has 
suffered irradiation driven mass loss. If we do not include 
illumination then P eq is longer, since at the beginning of 
RLOF mass transfer is proceeding on the thermal time-scale 
at a rather high mass accretion rate, and the neutron star 
spins-up to millisecond periods. After about 0.2 Mq have 
been lost, the mass transfer rate decreases and the neutron 
star enters a "propeller phase" during which the fast rotat- 
ing neutron star will slow down. To avoid irradiation-driven 
mass loss from this system the external irradiation flux must 
be less than critical flux ~ 1.6 x 10 10 ergs _1 cm -2 (Podsiad- 
lowski 1991). Let us estimate the value of the initial pulsar 
period Po to allow such a flux. Simple estimates show that 
Po must be longer than 3.85xl0~ 3 s, and the time elapsed 
from detachment of the Roche lobe is ~ 3.2 x 10 yrs. If Po > 
3.85x 10 _3 s then there is no irradiation driven mass loss and 
the hydrogen burning is stationary, prolonging the cooling 
time. Since both systems are in many parameters very simi- 
lar to each other, only the white dwarf cooling time-scale is 
very different and there must be something which leads to 
the different cooling histories for the helium white dwarfs in 
these systems. Different irradiation fluxes (F 1 a; (J0751+1807) 
?s 10 x Fr(J1012+5307)) may explain the discrepancy. 

We would like to point to an important issue neglected 
in many papers. Usually the comparison with observed sys- 
tem is made for one or several parameters (for example the 
cooling age of the white dwarf, T e f / and so on) but not with 
all available data. For binary systems we must build a self- 
consistent picture with the right orbital parameters (which 
may be a very hard task, as in the case of the two systems 
discussed here) , masses (if known) and cooling parameters 
of the white dwarf. 

The effect of irradiation on the secular evolution of bi- 
naries containing a millisecond pulsar has been discussed in 
a recent review paper by D'Antona (1996). She has pointed 
out that one needs to distinguish the two main effects of 
irradiation with following terms; i) illumination: structural 
changes in the whole stellar structure due to irradiation, 
which may cause the secondary to fill its Roche lobe again, 
ii) evaporation: direct driving of mass loss from the sec- 
ondary. Both cases give the same final result - the decrease 
of the mass of the hydrogen envelope. 

We propose that for PSR J0751+1807 the irradiation 
by the pulsar causes irradiation-driven mass loss from sec- 
ondary, resulting in exhaustion of the hydrogen envelope, 
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and as a consequence of this, a hydrogen shell burning 
stops. The irradiation-driven wind influences also the cool- 
ing phase, decreasing the cooling time-scale by a factor 1.7-4 
(depending on how efficient this process is) . 

For the PSR J1012+5307 system the irradiative flux 
may be too small to cause significant mass loss, and hydro- 
gen shell burning will be maintained for a long enough time 
to keep the white dwarf hot. 
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